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Oh ' In this paper, we deal with a class of backward doubly stochastic differential equa- 

^ cj tions (BDSDEs, in short) involving subdifferential operator of a convex function and 

driven by Teugels martingales associated with a Levy process. We show the exis- 
tence and uniqueness result by means of Yosida approximation. As an application, we 
give the existence of stochastic viscosity solution for a class of multivalued stochastic 
partial differential-integral equations (MSPIDEs, in short). 
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1 Introduction 

Backward stochastic differential equations (BSDEs, in short) related to a multivalued max- 
imal monotone operator denned by the subdifferential of a convex function have first been 
introduced by Gegout-Petit and Pardoux lfl4l . Further, Pardoux and Ra§canu Eol proved 
the existence and uniqueness of the solution of BSDEs, on a random (possibly infinite) time 
interval, involving a subdifferential operator in order to give a probabilistic interpretation 
for the viscosity solution of some parabolic and elliptic variational inequalities. Follow- 
ing, Ouknine [23], N'zi and Ouknine |[T9l l20l . Bahlali et al. [3]|4l discussed this type of 
BSDEs driven by a Brownian motion or the combination of a Brownian motion and an 
independent Poisson point process under the conditions of Lipschitz, locally Lipschitz or 
some monotone conditions on the coefficients. 
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Recently, a new class of BSDEs, named backward doubly stochastic differential equa- 
tions (BDSDEs, in short) involving a standard forward stochastic integral and a backward 
stochastic integral has been introduced by Pardoux and Peng [25] in order to give a prob- 
abilistic representation for a class of quasilinear stochastic partial differential equations 
(SPDEs, in short). Following it, Matoussi and Scheutzow IfTHll . Bally and Matoussi 0, 
Zhang and Zhao ll3T1l . Aman and Mrhardy[l] and Boufoussi et al. 1U|7) studied this kind 
of BDSDEs from different aspects. 

The main tool in the theory of BSDEs is the martingale representation theorem for a 
martingale which is adapted to the filtration of a Brownian motion or a Poisson point process 
(Pardoux and Peng (24), Tang and Li |[30l ). Recently, Nualart and Schoutens [21] gave 
a martingale representation theorem associated with a Levy process. This class of Levy 
processes includes Brownian motion, Poisson process, Gamma process, negative binomial 
process and Meixner process as special cases. Based on [21], they showed the existence 
and uniqueness of the solution for BSDEs driven by Teugels martingales associated with a 
Levy process in ||22|. These results were important from a pure mathematical point of view 
as well as from application point of view in the world of finance. Specifically, they could 
be used for the purpose of option pricing in a Levy market and related partial differential 
equation which provided an analogue of the famous Black-Scholes formula. Motivated by 
Il25l and IT221 . Ren et al. |[29l considered a class of BDSDEs driven by Teugels martingales 
and an independent Brownian motion, obtained the existence and uniqueness of solutions 
to these equations, which allowed to give a probabilistic interpretation for the solution to a 
class of stochastic partial differential-integral equations (SPDIEs, in short). Very recently, 
Ren and Fan |[28l derived the existence and uniqueness of the solution for BSDEs driven by 
a Levy process involving a subdifferential operator and gave a probabilistic interpretation 
for the solutions of a class of partial differential-integral inclusions (PDIIs, in short). 

Motivated by the above works, the first aim of this paper is to derive existence and 
uniqueness result to the following BDSDE involving subdifferential operator of a convex 
function and driven by Teugels martingales associated with a Levy process 

AY t + f(t , Y„ Z t ) At + g(t , Y t , Z t ) AB t G dq{Y t ) At + ZT= i zf dflf , < t < T, 

(1.1) 

where 3cp is a subdifferential operators. The integral with respect to {B,} is a backward 
Kunita-Ito integral (see Kunita lfl6ll ) and this one with respect to {Hj },->i is a standard 
forward Ito integral (see Gong lfl5l ). Our method is based on the Yosida approximation. 

On the other hand, since the pioneering paper due to Buckdahn and Ma [9],[10],[H ], 
the notion of stochastic viscosity solution has been intensely studied in the last ten year. 
Among others, we can cite the work of Boufoussi et al. 10, Q, Aman and Mrhardy ID, 
Aman and Ren [2] and Ren et al. 11271 . etc. In all these different works, authors have set 
existence results to stochastic viscosity solution of several types of SPDE. The tool is en- 
tirely probabilistic and used the connection between these SPDE and associated BDSDEs. 
Following this way, the second goal in this paper is to give stochastic viscosity solution for 
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multivalued stochastic partial differential-integral equations (MSPDIEs, in short) 



( $g(t,x) + Lu(t,x) +f(t,x,u(t,x), (u\ (t,x) )^ =1 ) +g(t,x,u(t,x))B t ) E 3(p(x), < K I, j 6 

(1.2) 

u(T,x) = wo(jc), x € R d , 
where £ is the second-order differential integral operator of the diffusion process given by 

L${t,x) = mil^a i (x) — (t,x) + -l^af{x)- 1 {t,x) 

1=1 OX < Z i=l 

+ f $(t,x + a(x)y)-y(t,x)-(yM,x),a(x)y)]v(dy), 

JR 

(1.3) 

and 

<t>lM = / (M*+<K*)y)-$M)/>*(y)v(dy), 



with a a M^-valued function, coefficient of SDE driven by the Levy process L, m\ = E(Li) 
and pk precise later . Our method is also fully probabilistic and uses connection between 
MSPDIE (0> and BDSDE (0> in Markovian framework. 

The paper is organized as follows. In Section 2, we introduce some preliminaries and 
notations. Section 3 is devoted to the existence and uniqueness result for BDSDEs involving 
subdifferential operator of a convex function and driven by a Levy process. Finally, in 
section 4 we derive a probabilistic representation (in stochastic viscosity sense) for the 
solution of a class of MSPDIEs via BDSDEs proposed in Section 3. 

2 Preliminaries and notations 

Let T > be a fixed terminal time and {B t : t € [0, 7 1 ]} be a standard M- valued Brownian 
motion defined on a complete probability space (O, f ,P). Let us also consider {L, : t 6 
[0, T]}, a R- valued Lev process corresponding to a standard Levy measure v, defined on a 
complete probability space (O', J ',P'), with the following characteristic function: 



E(e iuL ') = exp 



/awf - -K 2 u 2 t + t [ (e mx - 1 - £ieclj-| JC | <1 T.)v(dbc) 

2 JR 



where a € R,K > 0. Moreover, the Levy measure v satisfies the following conditions: 

1- / K (lAy 2 )v(d3;) <oo; 

2. Jj_ E E r c e^Mvfdy) < °°, for every £ > and for some X > 0; 

which provides that L, has moments of all orders, i.e. fj^ |x| ! v(dx) < °°, V/ > 2. 
Let consider the product space (H,5,P), defined by 

£2 = £2x£2'; 7=7®7'\ P = P<g>P'. 
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Further, random variables ^(co), CO € 12 and oo' € 12' can be considered as random 

variables on Q. via the following identifications: 

^(co,a)0 = ^co); £(oo,oy) = C(eo'). 

In this fact, the processes B and L are assumed independent. Next, denoting by 0\[ the 
totality of P-null sets of 7 , and for each t € [0, T] , we define 

where for any process {%}, = cr{ri r — rj s : s < r < t} and 7^ = 7q v Since {iF f ^-}r>o is 
decreasing and {!F t L }t>o is increasing, the objet {ft}t> is neither increasing nor decreasing. 
Thus it does not a filtration. 

We denote by (H^)/>i the linear combination of so-called Teugels martingale Y t as- 
sociated with the Levy process {L t : / 6 [0, T]}. More precisely 

H U=c i , i Y} i) +c U - l Yt 1) + --- + cuY} 1) , 

where for all i > 1, F f (0 = L ( W - £[4°] = L f (0 - r£[L^]. For each f € [0,T], L, {0 is a 
power-jump processes defined as follows: = Lt and Lj' = £o<s<f 

(AL S )' for / > 2. 

It was shown in Nualart and Schoutens ETTl that the coefficients c,-^ correspond to the 

orthonormalization of the polynomials (x) = Cijx l ~ l + Cjj-\x'~ 2 H hc, j with respect 

to the measure fj(dx) = x 2 v(dx) + K 2 8o(dx): 

/ q n (x)q m (x)fi(dx) = if n ^ m and / q 2 (x)fj(dx) = 1. 
Jm. Jw 

We set 

p k (x) =xq k -\(x). 

The martingales (ffW ),->i can be chosen to be pairwise strongly orthonormal martingales, 
i.e. [H^>,H^] = 0,i ^ j, and {[H^ ,H^] t — ?} f >o are uniformly integrable martingales 
with initial value and (H^,H^) { = 8 u t. 

Remark 2.1. The case of V = corresponds to the classic Brownian case and all non-zero 
degree polynomials qi{x) will vanish, giving H® = , i = 2, 3, • • ■ , i.e. all power jump 
processes of order strictly greater than one will be equal to zero. If v only has mass at 1 , we 
have the Poisson case; here also H, = 0, i = 2,3, • • • , i.e. all power jump processes will be 
the same, and equal to the original Poisson process. Both cases are degenerate in this Levy 
framework. 

Let us introduce the following appropriate spaces: 

• f = |x= (Xn)n>V, \\x\\p = (£" =1 |x n | 2 ) 1/2 < ooj. 

• rt 2 the subspace of the ^-measurable and M-valued processes (Y t ) te [ 0T ] such that 

rT 

\\Y\\ 2 H2 =E |F f | 2 d?<+oo. 
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S 2 the subspace of the R- valued, f t - measurable, right continuous left limited (rcll, in 
short) processes (if )?e[o,r] such that 



11711^=1 sup \Y t \ z <+oo. 

V0<(<:r / 

(P 2 (l 2 ) the space of jointly predictable processes (Z) fe [ r j taking values in I 2 such 
that 

Z||| 2(/2) = E /" Irtds^E /" |Z, w | 2 dj<oo. 

t/ ;_ j </ 



l2> 2 (/ 2 ) 

Now, we give the following assumptions: 



(HI) The coefficients / : [0, T] x Q x R x £ 2 -> R and g : [0, T] x £2 x R x f -)■ M satisfy, 
for all f € [0, 7] ,y G R and z € ^ 2 , 

(i) f(t,-,y,z) and g(t,-,y,z) are ^-measurable, 

(ii) /(-,0,0), g(-,0,0) G^ 2 ; 

(H2) There exists some constants C > and < a < 1 such that for every (t, go) G [0, T] x 

n,{yuzi),(y2,z 2 )eMxf 

\f(t,(o,y u zi)-f(t,(0,y 2 ,z. 2 )\ 2 <C(\ yi -y 2 \ 2 + \\zi -z 2 \\%) , 

\g{t,<ji,y u zi)-g{t } <Si,y 2 ,z 2 )\ 2 <C\yi -y 2 \ 2 + a\\z\ -z 2 \\ 2 2 ; 

(H3) Let cp : R — » (— °°, +°°] be a proper lower semi continuous convex function satisfying 
q>(y) > 9(0) = 0; 

(H4) The terminal value £ G L 2 (£l, jTr,P) satisfies 

E(|^| 2 + cp©)<oo. 

Define: 

Dom(cp) = {«Gl: q>(«) < +00} , 

3cp(w) = {u* GR : (w*,v-w) + cp(w) < cp(v), for all v G R}, 

Dom(3<p)) = {u G R : 9cp(w) 7^ 0}, 

Gr(3q>)) = {(«,«*) G R 2 : w G Dom(3cp),w* G d<p(u)}. 

Now, we introduce a multi-valued maximal monotone operator on R defined by the subdif- 
ferential of the above function (p. 
For all x G R, define 

cp E (jc) =minQ|x-.y| 2 + £(p(.y)^ =^\x-J £ (x)\ 2 + ey(J e x), 

where J e (x) = (I + e3cp) — 1 (jc) is called the resolvent of the monotone operator A = Bcp. 
Then, we have the following Proposition which appeared in Brezis [8 ]. 
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Proposition 2.2. (1) (p e : R — > R is a convex function with Lipschitz continuous deriva- 
tives; 

(2) Vx G R, ±D<Pe(*) = 53<fc(*) = \{x-J z (x)) G acp(/ E (x)); 

(3) Vx,yeM,|7 E (x)-7 E (3;)| < \x-y\; 

(4) Vx € R,0 < (p e (x) < (Dcp E (x),x) ; 

(5) Vx,y G R W £,5 > 0, (±Dcp E (x) - iDcp 8 (3;),x-3;> >-(* + $) |Z><p E (*)ll»<PsO0l- 

We first give the definition of BDSDEs involving subdifferential operator of a convex 
function and driven by Levy process. 

Definition 2.3. We call solution of BDSDE (^,/,g,cp) atripleof (Y,U,Z) of progressively 
measurable processes such that 

1. (F,Z) eS 2 x 2> 2 (/ 2 ), U G ^ 2 ; 

2. (F t ,Z7 t ) e3<p, dP(g)d?-a.e. on [0,T]; 

/■T pT pT 00 /*T 

3. Y t + y U s ds = 1~ + j f{s, Y s , Z,) ds + J g(s, Y s , Z,) dB s - £ y zj° dffj , < f < 

r. ' ' 

3 Existence and uniqueness result for BDSDE driven by Levy 
process 

The first result of the paper is the following theorem: 

Theorem 3.1. Assume the assumptions o/(Hl)-(H4) hold. Then, the BDSDE (^,/,g,cp) 
has a unique solution. 

For the prove of this theorem, let us consider the following BDSDEs: 

y f E + ^ r Dcp E (F/)ds = $ + jT 7 /(s,F s E ,Z E )ds + J*g(s,Y?,Z*)dB, 

-£ / z^dffj , o<? <r, (3.1) 

(=i 

where cp E is the Yosida approximation of the operator A = 3cp. Since ^D(p E (T s E ) is Lipschitz 
continuous, it is known from a recent result of Ren et al. [29], that Eq. (13.11 ) has a unique 
solution (T E ,Z E ) G S 2 x !P 2 (/ 2 ). 

Setting £/f = ±Dq> e (y, E ), < t < T, our aim is to prove that the sequence (F E ,£/ E ,Z E ) 
converges to a sequence (Y,U,Z) which is the desired solution of the BDSDEs. 

In the sequel, C > is a constant which can change its value from line to line. Firstly, 
we give a prior estimates on the solution. 
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Lemma 3.2. Assume the assumptions of (H1)-(H4) hold. Then, there exists a constant 
C\ > such that for all £ > 

E( sup |F f e | 2 + f T \\Z*\\ 2 f2 ds) <d. 

Vo<?<T ^0 / 

Proof. Applying the Ito formula to \Y t e \ 2 yields that 

\Y t E \ 2 + lf t T Y*D<? e (Y*)ds = n 2 + 2f\f-f(sJ!^ s )ds + 2f\^g{sjf^)dB s 

+ 1 \ 8 {sJ!^)\ 2 ds-^J t |Zf«| 2 d[H« H«] s 

-2£ / Y^dHP. (3.2) 



i=i JT 

Noting that the fact F s E Dcp e (y/) > and taking expectation on the both sides, we obtain 

rT rT 

E|7 f e | 2 + E^ \\Z*\\%ds < E|^| 2 + 2E^ F 4 E -/(s,F/,Z s e )ds 

+Ej t T \g(s,Y s E ,Z*)\ 2 ds. (3.3) 

Using the elementary inequality lab < $ 2 a 2 + & for all a,b>0 and (H2), we get 

2yf(s,y,z) = 2y(f(s,y,z)-f(s,0,0)) + 2yf(s,0,0) 



< 



1 

M 



\y\ 2 +MC\y\ 2 +MC\\z\\% + \y\ 2 + |/M,0)f 



< ( 1 + ^+MCj |v|-H/|,.0.0)| J +A/r||;;| l ,. : . 



\g(s,y,z)\ 2 = \g(s,y, z )-g(s,0,0)+g(s,0,0)\ 2 

< + \g(s,y,z)-g(s,0,0)\ 2 + (1 + P)k(s,0,0)| 2 

< (i + cm 2 + (i + P)kM,o)| 2 + « (i + 1) II*- 

Choosing M = f3 = it follows from (O that 

E|F, £ | 2 + ^E| r ||Zf|| 2 2 d, 

n 2 + f t \Y!\ 2 ds + \f(s,0,0)\ 2 ds + J o \g(s,0,0)\ 2 ds 
Gronwall inequality and Bulkholder-Davis-Gundy inequality show the desired result. □ 
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Lemma 3.3. Assume the assumptions o/(Hl)-(H4) hold. Then, there exists a constant 
C2 > such that 

(OE^ T Q|D(p E (F s e )|^) ds<C 2 ; 

(ii) E(p(7 e (F ; E ))<C 2 ; 

(n7)E|F f E -7 e (y f E )| 2 <s 2 C 2 . 

Proof, (i) Given an equidistant partition of interval [0, T] such that = to < ?i < t 2 < ■ ■ ■ < 
t„ = T and ti + \ — fj = K the subdifferential inequality shows 

cp E (F, E +1 ) > cp E (F E ) + (Y^ - F, E )Dcp E (F f E ). 

From (13 -lb . we obtain 

+ I jT 1 Dcp E (F r E )Dcp E (F 4 E ) d, < cp E (F r E +i ) + £ +1 D(p E (F f E )/(,, F/,Z E ) dj 

Jti 

-2£/' D(p e (F ; E )(Z E )^d//, W . 
Summing up the above formula over i and letting n — > °°, we obtain 

<?e(Y t t) + y o T \D<p e (Y*-)\ 2 ds 

< cp e ©+ [ Dq>z(Y s e -)f(s,Y s e ,Z*)ds+ f D<p e (Y*-)g( S ,Y s e ,Z%)dB s 
Jo Jo 

-2£ / Dcp e (F E )(Z E )°' ) d//i i) 

Hi JO 



Taking expectation on the both sides, we get 



Ecp E (F f E ) + -E / T |Dcp E (F 4 E )[ 2 ds<Ecp e (^)+E f D^pfiffaYf^ds. (3.4) 

8 JO JO 



For 



D<p E (y)f( S ,y,z) < ^\Dq> E (y)\ 2 + ^\f(s,y,z)\ 2 



< Yz I 2 + £ ( l/fr* ^) - /(■*> o,o)| 2 + 1/(5, 0, 0) | 2 ) 

< ^|D(p E (3;)| 2 + sCb| 2 + sC||z|| 2 2 +s|/(^0,0)| 2 , 



the fact that cp e (F r E ) > and cp E (^) < ecp(^), we obtain 

1 T / T T 

-E[ |Dcp E (F/)| 2 d.<CE((p©+ f |/(.,0,0)| 2 d. + r sup \Y t £ \ 2 +f \\Z?\\ 2 2 dt 
2e Jt V Jo o<t<T Jo 

Lemma l3T2l shows the desired result, 
(ii) From (13.41 ). we obtain 

Ecp e (F, E ) < elcp© + \D^{Yf)\ 2 ds + €&J* \f(s,Y?,Z*)\ 2 ds. 

Using (p(/ e (F f E )) < ±(p e (F f E ) and (i), we obtain (ii). 

The last part of the Lemma simply follows from the fact that 

|x-/ e (x)| = 2(p e (x) -2£cp(7 e (x)). 

□ 

In what follows, we aim to show that (F E ,Z E ) is a Cauchy sequence in S 2 x $ (I )■ 

Lemma 3.4. Assume the assumptions o/(Hl)— (H4) hold. Then, there exists a constant C3 
such that for all e, 8 > 

e( sup |F f E -F f 8 | 2 + f T \\Zf-Zf\\% dt] <C 3 (£ + 5). 
\o<f<r Jo / 

Proof. Applying the Ito formula to |F, E — F 8 | 2 yields that 

| F e_ F 8|2 = _2^ r (F/-F/)QD( Pe (F/)d J -iD(p8(j; 5 ))d, 
+2 j* (F S E - F s 8 ) (f(s, F/,Z e ) - f(s, F/, Zf ) ) ds 
+2^ r (F s e -F s 8 )fe(,,F/,Z E )-g(,,F s s ,zf))dB 1 
+ £ \g(s^^)- g (sj?,zf)\ 2 ds- £ Itf-Z^ds 
-2£/ (Y s E -Y s 5 )(Z^' (l) -Z^)dHP. (3.5) 

Taking expectation, we obtain 

E|F f E - F r 8 | 2 + EJ ||Z E - Zf|| 2 2 ds 

= -2Ejr r (l?-r/)Qz^(F/)dr-iD9 S (y, 8 )) ds 
+2E^ r (F E -F $ 8 )(/(,,F 4 E ,Z E )-/(,,F 8 ,zf))d, 
+E jT |^,F/,Z e ) -g(s,F 8 ,Z 8 )| 2 ds. (3.6) 
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Using the elementary inequality lab < $ 2 a 2 + & for all a,b > and (H2), we get 
{Yf-Y*){f(s,Y*,Z*)-f{ S jf,Z*)) 



2C i y e v 5|2 , 1_a i v e v 8|2 , 1 - °W .yO N 2 

1-oc 1 s 



< Itf " + -r 2 I*? - i?r + -r 2 113 " tine 



and 

|^,F/,Zf)-g(,,y8 )Z f)|2< C | y/ _ F /|2 + a || zE _ Z S||2 2 _ 

Noting (5) of Propositio ri2.2l we obtain 

nY ?-Y?\ 2 + l -^Ef t T \\Zt-Z%ds 

< CeE [ T \Y?-Y*\ 2 ds 



+2Q + i)E| r |D(p E (y s E )||D(p 8 (y s 5 ))|d,. (3.7) 

Lemma |331 shows that 

2 {l + l) E f t T l D( P E ( F /)H D( Psft S ))|ds < (e + 8)C. 

So, we can obtain 

E|r e -r f 8 | 2 + Ey ||Z E -Z 4 8 || 2 2 ds<CEy |F/-F 8 | 2 ds + C(£ + 8). 



The Gronwall inequality shows that 

E|F E -F, 8 | 2 + E f T 

0<t<T 

The Bulkholder-Davis-Gundy inequality shows the desired result. □ 



sup E|y f e -y f s | 2 + E / \\Z*-Z*\\%ds<C(e + S). 



Proof of Theorem 13.1 1 

Existence. Lemma [3~4l shows that (F E ,Z E ) is a Cauchy sequence in S 2 x P 2 (l 2 ). Denoting 
its limit by (F,Z), then it follows from Lemma 1331 (F,Z) € S 2 x V 2 (l 2 ). For each £ > 0, 
define Uf = ±Dcp £ (F, E ) and fjf = ft Ufds. Therefore d3TH) and Lemma El yield that for all 
£, 8 > 

E ( sup 10* - 0?\ 2 J < CE ( sup \Y t e - F 8 | 2 + [ T \\Zf - zf\\% At] , 

\0<t<T J \0<t<T JO J 

which shows that (f/ e ) is a Cauchy sequence. Hence, there exists a measurable process U t 
such that 

limE( sup \U?-U,\ 2 ) =0. 
Vo<f<r / 
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Furthermore, Lemma 1331 (?) shows that 



supE [ T \U?\ 2 dt = supE [ T ( -\D<p e (Y t £ )\) dt 

£ JO £ JO \£ J 



which shows that Uf is bounded in the space L 2 (D.,H l [0, T]), and (U E ) e converges weakly 
to a limit in that space and the limit is necessarily U. In particular, U is absolutely continu- 
ous. So, there exists a measurable process (U t )o<t<T G ^ 2 such that U t = /gL^ds. 

Next, we show that (Y t ,U t ) G 3<p, dP(g> df-a.e. on [0,T]. Moreover, with the help of 
Lemma 5.8 in QU, and for all < a < b < T, V G ^ 2 ([a,Z?]), we obtain 



f b Uf{V t - Y t E ) dt j h U ; (V t - Y t ) dt, as £ -> 



in probability. In particular we have 



[ Uf(J e (Y t e ) - Y t e ) dt 0, as £ -> 0. 
which together with Proposition l2.2l provides that £/, E G 9(p(7 E (F f E )) and 
[ b U?(V t -J e (Y t E )dt+ f\(J £ (Y t e ))dt< f\(V t )dt. 

J a J a J a 

Taking the liminf in probability in the above inequality, we obtain 

f b U t (V ! -Y ! )dt+ [\(Y t )dt< [\{V t )dt. 

Ja Ja Ja 

Since a,b and the process V are arbitrary, this shows that 

Ut(V t -Y t ) + y(Y t ) < q>(V J ),dP<g>d/-a.e. on [0,r]. 



Taking limit on the both sides of (13.11) . we obtain the existence of the solution. 

Uniqueness. Let (Y t ,U t ,Z t )o< t <T and (Y/ ,U{ ,Z' t )o<t<T be two solutions of BDSDEs 
associated with (^,/,g,cp). Define 

(AY t ,AU t ,AZ,)o< t < T = (Yt-Y^Ut-U^Zt-ZDo^r. 

Applying the Ito formula to \AY t \ 2 shows that 

rT rT 

E|AF,| 2 + 2E J AU S AY S ds + ^J t 11^ III ds 

= 2Ej\Y s [f( S ,Y s ,Z s ) -f(sX,Z' s )]ds 

\g(s,Y s ,Z s )-g(sXX)\ 2 &- (3-8) 

Since 3cp is monotone, we obtain 

AU S AY S > 0, dP®df-a.e. 
Further, as the same procedure as Lemma |3~4l we obtain 

rT rT I r-T 

E|Ay,| 2 + E^ ||AZ ( || 2 2 ds<CEy |Ay 4 | 2 ds + -Ey ||AZ 4 .|| 2 2 ds. 
The Gronwall inequality shows the uniqueness of the solution. 
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4 Stochastic viscosity solutions of multivalued SPDIEs 

In this section, we derive the existence of the stochastic viscosity solution of a class of mul- 
tivalued SPDIE (11.21) via BDSDE with subdifferential operator and driven by Levy process 
studied in the previous section. 

4.1 Notion of stochastic viscosity solution of multivalued SPDIEs 

Let us recall F B = {^^-}o<f<r be the filtration generated by B. The objet WCq T denotes all 
the F B -stopping times x such < x < T , a.s. and M 3 is the set of all almost surely finite 
F B -stopping times. For generic Euclidean spaces E and E\, we state those spaces: 

1. The symbol C k,n ([0,T] x E;E\) stands for the space of all E\ -valued functions de- 
fined on [0, T] x E which are /c-times continuously differentiable in t and «-times 
continuously differentiable in x, and C^"([0,T] xE;E\) denotes the subspace of 
C k, "([0,T] x E;E\) in which all functions have uniformly bounded partial deriva- 
tives. 

2. For any sub-G-field g C rf, c k ' n (g,[0,T] XE;E X ) (resp. c£' n (g , [0, T] xE;E{)) 
denotes the space of all C k ' n ([0, T] x E;E\) (resp. C^"{[0, T] x E;E\) -valued random 
variable that are g ® S ([0, T] x £')-measurable; 

3. C k ' n {F B ,[0,T] xE;E{) (resp.cf "(F 5 , [0,T] xE;Ei)) is the space of all random fields 
cp G C hn (j T , [0,T] x E;E X (resp. C Kn (f T , [0,T] x E;E y ), such that for fixed x G E 
and t G [0,T], the mapping go — > a(t,(0,x) is F B -progressively measurable. 

4. For any sub-G-field g C f 3 and a real number p>0,L p (g;E) denotes the set of all 
^-valued, (^-measurable random variable £, such that < °°. 

Furthermore, regardless of the dimension, we denote by (•, •) and | • | the inner product and 

norm in E and E\ , respectively. For (t,x,y) G [0, T] x W 1 x R, we denote D x = (3775 •> 5^7)' ^xx = 

0xix-)ij=i> Dy = j-, A = The meaning of D^. and D >v is then self-explanatory. 
The coefficients 

ilx[0,T]xR 1 'xlxf 2 -fK 
OxfOjrJxK^xR-)! 
R (/ ->• R rf 
M rf — >■ R, 

satisfying assumptions: 

f <tf(l + |*| + |y| + ||z||), 

(H5) 

[ \u (x)\ + \<?{u (x))\ <K(l + \x\). 

( \\a(x)-a(x')\\ <K\x-x'\, 
(H6) \ 

[ \f(t,x,y,z)-f{t,x,yj)\ < K(\y-y'\ + \\z-z!\\p). 
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(H7) The function g G C°' 2 ' 3 ([0,r] xl^x R;R). 

The definition of stochastic viscosity solution to MSPDIE dl.ll ) use the stochastic sub- 
and super-jets introduced by Buckdahn and Ma f9). Let us recall the following needed 
definitions. 

Definition 4.1. Let x G Mq T , and £ G f x . We say that a sequence of random variables 
(Xk,^k) is a (x, ^-approximating sequence if for all k, (x*,^) G x L 2 (j x ,R rf ) such that 

(i) ^ — ^ £, in probability; 

(ii) either x^ f X a.s., and x^ < X on the set {x > 0}; or Xt 4- X a.s., and x k > x on the set 

{x<r}. 

Definition 4.2. Let (x,£) G fW B r x L 2 ( J t B ;M d ) and u G C (F B , [0, T] x M d ) . We denote by 

Jg ,2 ' + u(i, Q the stochastic g-superjet of u at (x, the set of R x R d x 5 («)-valued and 5 X B - 
measurable random vector (a,p,X) (S (d) is the set of all symmetric d xd matrix) which is 
such that for all (x, ^-approximating sequence (x^,^), we have 

u{x k £ k ) < u{T,Q+a(T k -T) + b(B Zk -B z ) + ^(B Zk -B z ) 2 + (p£ k -ty 

+(q&-z.) -B x ) + l(x^ k -^ k -t,) 

+0 (| Tft _ T |) +0 (|k_£[2). (4.1) 
The ^"/-measurable random vector (b,c,q) taking values in R x R d x R d is defined by 

< 

k q = d x g(%&,u('C,Z ) )) + d u g(i,£ ) ,u('Z,t,))p. 

Similarly, jg ,2 ~u(T,Z,) denotes the set of all stochastic g-subjet of u at (x,^) if the inequality 
in (14.11 ) is reversed. 

RemarkA3. Let us note that 3cp(y) = [(pj(y),<pj.(y)], for every y G Dom(<p), where tpj(y) and 
(pj.(j) denote the left and right derivatives of (p. 

In order to simplify notation in the definition of the notion of stochastic viscosity solu- 
tion of multivalued SPDIEs, we set 

V f (<z£,a,p,X) = - a -^Tmce(aa*^)X)- mi (p,a^))-U (Xa(^),a(^))v 2 v(dy) 

-f(x^M^)J & (pM^y)Pk(yMdy)Y 

Definition 4.4. (1) A random field u G C (F B , [0, T] x R rf ) which satisfies u (T,x) = u (x), 
for all x G R rf , is called a stochastic viscosity subsolution of MSPDIE (11.11 ) if 

«(t,$) G Dom(cp), V (x,£) G Mlj x L 2 (V T B ;M d ) , P-a.s., 
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and at any (%&) G x L 2 (f x B ;R d ), for any (a,p,X) G j/^ 2 ' + w(x,^), it hold P-a.s. 

Vf(xX,a,p,X) + (p[(u(^)--(gd u g)(^,u(^)) <0; (4.2) 

(2) A random field u E C (F B , [0, T] x M d ) which satisfies that w (r,x) = u Q (x), for all x G 
R' 5 ', is called a stochastic viscosity supersolution of MSPDIE d 1 -2b if 

«(t,$) G Dom((p), V (x,£) G 3W B r x L 2 (V T B ;R rf ) , P-a.s., 
and at any (x,£) G fW B r x L 2 (jF T B ;R d ), for any (a,p,X) G y g 1,2 "m(x,^), it hold P-a.s. 

V(T^,a,p,X) + cp;( M (x,^) - l -(gd u g)(x^u(^)) > 0; (4.3) 

(3) A random field u G C (F B , [0, T] x R rf ) is called a stochastic viscosity solution of MSPDIE 
(11.11) if it is both a stochastic viscosity subsolution and a stochastic viscosity supersolution. 

Remark 4.5. Observe that if / is deterministic and g = 0, Definition 14.41 becomes the gen- 
eralization of the definition of (deterministic) viscosity solution of MPDIE given by N'zi 
and Ouknine in ll20ll . 

To end this section, we state the notion of random viscosity solution which will be a 
bridge link to the stochastic viscosity solution and its deterministic counterpart. 

Definition 4.6. A random field u G C(F B , [0, T] x W 1 ) is called an co-wise viscosity solution 
if for P-almost all go G £1, w(co, •, •) is a (deterministic) viscosity solution of MSPDIE (1 1 -2b - 



4.2 Doss-Sussmann transformation 

In this section, using the Doss Sussmann transformation, our aim is to convert a multi- 
valued SPDIE to a PDIE with random coefficients so that the stochastic viscosity solution 
can be studied co-wisely. We first establish the link between the g-super or sub jet of u the 
solution to multi-valued SPDIE and super or sub jet of v solution to the converter PDIE with 
random coefficients. For instance let us consider the stochastic flow r\ G C(F B , [0, rjxl^x 
R), unique solution of the following stochastic differential equation in the Stratonovich 
sense: 

r\(t,x,y) = y + J (g(s,x,r\{s,x,y)),odB s ), (4.4) 

where (14.41) should be viewed as going from T to t (i.e y should be understood as the initial 
value). Under the assumption (HI), the mapping y i-> r\(t,x,y) defines a diffeomorphism 
for all (t,x), P-a.s. such that its y-inverse e(t,x,y) is the solution to the following first-order 
SPDE: 

^(t>x,y) = y- J (D y e(s,x,y),g(s,x,r\(s,x,y))odB s ). 
We refer the reader to their paper ifTUl for a lucid discussion on this topic. We have 
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Proposition 4.7. Assume that the assumptions (H1)-(H7) hold. If for (x,^) G M B T x 
L 2 {j B ;R d ),u£C (F B ,[0,r] xR d ) and(a u ,X u ,p u ) belongs to J g' 2 ' + u(l,Z,), then (a v ,X v ,p v ) 
belongs to Jq 2,+ v(x,%), with v(-, •) = £(•, ■,«(•, •)) and 

a v = D > .e(t,£,m( t ^)K 
p v = D v £(x, 5, w(x, + D x e(x, u(t, ) 

X v = D } .£(x,^ M (x^))X H + 2D^£(x,^ M (x^)) j p* 
+ D xc e(x, ii (t, £)) + D v:y £(x, 4, • 

Conversely, if for (x,£) G SK,f r x L 2 {T*\R d ), v G C (F B , [0,T] x M d ) and (a v ,X v ,/? v ) G 
J7 ! ' + v(x,^), */ien (ou,X u ,p u ) G V ,+ m(x,£) wzY/z w( v ) = r|(-, •, v(-, •)) and 

a u =D y X\(x&,v(%&))a v 

p u =D y x\{%&,v{%&))p v +D x ^(%£,v(%&)) 

X u = D y r\ (x, v(x, 4) )X V + 2D xv ti (x, ^, v(x, %)) P ; 
+ D„r| (x, v(x, £)) + D vv ri (x, v(%£)) PvP * . 

However, contrary to classical SPDE, the resulting PDIE from MSPDIE (11.21 ) due to 
Doss-Sussman transformation is not necessarily MPDIE studied by N'zi and Ouknine (see 
GUI ). Therefore, we need the following version of viscosity solution for resulting PDIE 
obtain by Doss-Sussman transformation. 

Corollary 4.8. Assume that the assumptions (H1)-(H7) hold. Let us define and consider 
(%£) G Mlj x L 2 (F x B ;R d ), uec (F B , [0,T] x R d ). 

(1) for (a u ,X u ,p u ) G Jg' 2,+ u(T,t,)> u satisfies A4.2\l if and only ifv(-, •) = £(•, •,«(•, •)) satis- 
fies 



V^(x^,a v ,p v ,X v ) + 



D,ti(x,^v(x,4)) 



<0; 



(4.5) 



(2) for (a u ,X u ,p u ) G Jg' 2, «(x,^), u satisfies j[4. 3\) if and only ifv(-,-) = &(•,•,«(•,•)) satis- 
fies 



Vf(x,^,a v ,p v ,X v ) + 



cp;(Tl(x,^,v(x^)) 



where (a v ,p v ,X v ) is defined by Proposition \4. 7\ and 



f(t t x,y,(&)*>i] 



1 



>0; 



/ (t,x,T}(t,x,y),D y T}(t,x,y)Q k +T] l k (t,x,y) 



(4.6) 



Dyi\(t,x,y) L- 

- j Tl (*>*»?))+ M (*>*o0 + X(a* (x)D^,ri (f ,x,y), a(x) Pv [ 

+-XDyyT\(t,x,y)\c{x) Pv \ 2 



withQ k = f R (p v ,c(x)u)pk(u)v(du) and\= 1 + J R u 2 v(du) 
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Proof. Let G 5W B r xL 2 (^ T B ;lR rf ) be given and (a u ,p u ,X u ) G j7g' 2 ' + w(T,^). Weassume 
that m is a stochastic subsolution of MSPDIE (1 1 -2b . i.e. 



11(1;,$) G Dom(cp), V (x,t,) £M { ? T xL 2 [ff;R d ) , P-a.s., 

such that 

V / (T 1 tfl ) A^) + <P!("(x,5)-^(s3 1( g)(x ) § 1 «(T^)) <0, P-a.s. 

In view of Proposition 14.71 and since D y T[(t,x,y) > 0, for all (t,x,y), we obtain by little 
calculation 



q/ f (Tl(T,tv(T^))) 

D V T1(T,^V(X^)) 



Vy(x,^,a v ,j7 v ,A v )+ - & g ^ <0. 



The converse part of (1) can be proved similarly. In the same manner one can show the 
second assertion (2). □ 



5 Probabilistic representation result for stochastic viscosity so- 
lution to MSPDIEs 

In this section, we aim to show that the solution of multivalued BDSDE with jump gives 
the viscosity solution of a semi-linear MSPDIE in the Markovian case. 

5.1 A class of reflected diffusion process 

We now introduce a class of diffusion process. Let a : M. d — > M. d be a uniformly bounded 
function satisfying the uniform Lipschitz condition with some constant C > 0, for all x,y £ 

R d : 



\o(x) —o(y)\ < C\x — y\ 



(5.1) 



For each (t,x) G [0, T] x M. d , from [ 17 ] and reference therein, let {Xl' x ,s G [t, T] } be a unique 
pair of progressively measurable process, which is a solution to the following stochastic 
differential equation: 



X t / = x + J c{X' r L x )dL r . 
Furthermore, we have the following proposition. 



(5.2) 



Proposition 5.1. There exists a constant C > such that for all0<t <t' <T and x, x' G R d , 
such that 



E 



sup 




4" 


_0<s<T 







<C(\t'-t\ 2 + 



J\4 



X — X 



) 
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5.2 Existence of viscosity solution for MSPDIEs 

Fix T > and for all (t,x) 6 [0, T] x W 1 , let X'/ , s G [f, T] denote the solution of the SDE 
(15.21) . And we suppose now that the data (£,,f,g) of the multi- valued BDSDE with jump 
take the form 

£ = «o(Xr X )> 
/(j,y,z) = /W,y,z), 
g(s,y) = f(s,Xl*,y). 

And we give the following assumptions: 

We assume that u Q G C(R d ; R) , / G C( [0, T] x R J x R x £ 2 ; R) and g G C( [0, T] x R J x R; M) 
such that assumptions (H1)-(H7) hold. It follows from the results of the Section 3 that, for 
all (t,x) G [0,r] x W 1 , there exists a unique triplet (Y l ' x ,Z t ' x ,U t ' x ) for the solution of the 
following 

(1) (Y'^,Ul' x ) G 3cp, dP® ds,-a.e. on 

(2) F/>* + [ T U t /dr = u Q (X t /)+ f f{r,X t r X X/,Z t /)Ar+ f g^X'/X/W,. 

J s J s J s 

CO -J* 

-£ / {z*^m^,t<s<T. 

We extend processes F f A , Z** t/*** on [0, T] by putting = Y,'' x , Z 1 / = 0, iff = 0, 5 G 
[0,*]. 

We have this result whose proof is similar to that of Theorem 2. 1 appear in [ 25 ] 

Proposition 5.2. Let the ordered triplet (F/'*, Zj'*) be the unique solution of the multi- 
valued BDSDE (f5T2b . Then, for (s,t,x) G [0,T] x [0,T] x W 1 , the random field (s,t,x) ^ 
K'(Ys' x ) is a.s. continuous (Y t,x has jump), where E' is the expectation with respect to 
introduced at page 3. 

We are ready now to derive our main result in this section. 

Theorem 5.3. Assume the assumptions (H1)-(H7) be satisfied. Then, the function u(t,x) 
defined by u(t,x) = F/'* is a stochastic viscosity solution of MSPDIE (11.21) 

Proof. The fact that u(t,x) = Y t * = W(Y t "*), does not depend on go', it follows from Propo- 
sition IS2] that u G C{f B , [0,T] x R d ). Next, for all (x,£) G M B {0,T) x L 2 (> B ,R rf ), 

9(«(x(Q))^(a)))) = y(Y^ m ) <oo, P- a .s., 

which implies that m(t,£) G Dom(cp) P-a.s.Thus it remains to show that u is the stochastic 
viscosity solution to MSPDIE (11.21) . In other word, using Corollary 14. 81 it suffices to prove 
that v(t,x) = e(t,x,u(t,x)) satisfies (14.51 ) and (14.61) . In this fact, for each (t,x) G [0,T] x 
e, 8 > 0, let {(Ys^X^), 0<s<T} denote the solution of the following BDSDE: 

7 ;^8 + l^ r D(p5(y; .,,8 )dr = UQ (x t /) + ^{rX^Yf^K^Ar 

+ f\{rX- X X^)&B r 

CO „f 

-£ / (Z^ 8 )^d^ r W . (5.3) 



17 



Setting Y t '' x ' = u?(t,x), it is shown by Theorem 3.6 in [2], that the function v^(t,x) = 
e(t,x,u 5 (t,x)) is an co-wise viscosity solution to this MSPDIE: 

(i) (^M- [Lv 8 (t,x)+Mt,x,v 8 (t,x),o*(x)Vv & (t,x))\ j =0, (t,x)e[0,T]xR d , 



(5.4) 



(ii) v(T,x) = uo(x), x G M. d , 
where 



fb(t,x,y,z) =f(t,x,y,z) 



^Dy 5 (r\(t,x,y)) 
D } .ri(f,^,v) 



Moreover, letting Lemma 3.5 provide, along a subsequence that 

\v s (t,x) —v(t,x)\ -> 0, a.s., as 8 — ► 0, 



(5.5) 



for all (t,x) € [0,r] x R rf . 

On the other hand, for all (x,£) G Stf B (0,r) x L 2 (^ B ,M J ) and oo G £2 be fixed, let 
consider (a v ,p v ,X v ) G j7 1,2 ' + (v(x(co),i;(co))). Thus, since v 8 is an co-wise viscosity solution 
to the MSPDIE (15 -4b . and by Crandall- Ishii-Lions in [ 12], there exist sequences 

f 8„(gd)\0, 

(x„(co),^(co))G[0,r]xM rf , 

I «X>*v) e io' 2,+ ( v8n (^(a)),^(G)))) 
satisfying 

(x„((o),^(co),<,^,X v ",v 5 "(T„(a)),^((o))) (x(oo),^((o),a v , j p v ,^v,v(x((o),^((o))) 



such that for (x„(co),£„(co)) G [0,T] x 



V~ (aj) (x„ (co) , (co) , < , Z v n ,/>*) + 



« v » ^ i ZH Pai.( T l( T »( a> )'^( a) )' v5, (^( (D )'&»( (ll> )))) 



D V T1 (X„ (CO) , \ n (CO) , V 8 " (x„ (CO) , (q) ) ) 



<0. (5.6) 



In order to simplify the notation, we remove the dependence of CO. Let y G Dom(cp) such 
that y < «(x,^) = r|(x,^, v(x,^)). The ucp convergence of v 8 " to v implies that there exists 
«o > such that V« > «o, y < r)(x„,^„,v 8 "(x n ,^,j)). Therefore, inequality (15.61 ) yields 

(ti(x ;i ,^,v 8 "(x ;i ,^)) -y) V k (x n MX, P :) 

cp(v)-cp(7 8 „( 1 l(x^,v 8 «(x^))))l 



< 



D y X\{x n Xn,V h "{X n ^n))' 

Taking the limit in this last inequality, we get for all y < T|(x,^, v(x,^)) 

«P(Tl^>^v(x,^)))-9(y) 1 
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which implies that 



v . e v v , (p;(-n(T,^v(T,^))) 

^§,0,,^) + -^^^ <0, 

and derives that v satisfies (14.5I ). Hence, according to Corollary 14.81 u is a stochastic viscos- 
ity subsolution of MSPDIE dl .21) . By similar arguments, one can prove that w is a stochastic 
viscosity supersolution of MSPDIE d 1 -2b - This completes the proof. □ 
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